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PREFACE 


The  Electromagnetic  Compatibility  Analysis  Center  (ECAC)  is  a Department  of 
Defense  facility,  established  to  provide  advice  and  assistance  on  electromagnetic 
compatibility  matters  to  the  Secretary  of  Defense,  the  Joint  Chiefs  of  Staff,  the  military 
departments  and  other  DoD  components.  The  Center,  located  at  North  Severn,  Annapolis, 
Maryland  21402,  is  under  executive  control  of  the  As»stant  Secretary  of  Defense  for 
Communication,  Command,  Control,  and  Intelligence  and  the  Chairman,  Joint  Chiefs  of 
Staff,  or  their  designees,  who  jointly  provide  policy  guidance,  assign  projacts,  and  establish 
priorities.  ECAC  functions  under  the  direction  of  the  Secretary  of  the  Air  Force  and  the 
management  and  technical  direction  of  the  Center  are  provided  by  military  and  civil  service 
personnel.  The  technical  operations  function  is  provided  through  an  Air  Force  sponsored 
contract  with  the  IIT  Research  Institute  (IITRI). 

This  report  was  prepared  for  the  Systems  Research  and  Development  Service  of  the 
Federal  Aviation  Administration  in  accordance  with  Interagency  Agreement 
DOT-FA70lfVAI-175,  as  part  of  AF  Project  649E  under -Contract  F-19628-78-C-0006,  by  the 
staff  of  the  IIT  Research  institute  at  the  Department  of  Defense  Electromagnetic 
Compatibility  Analysis  Center. 

To  the  extent  possible,  all  abbreviations  and  symbols  used  in  this  report  are  taken  from 
American  Standard  Y10.19  (1967)  “Units  Used  in  Electrical  Science  and  Electrical 
Engineering"  issued  by  the  USA  Standards  institute. 
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FEDERAL  AVIATION  ADMINISTRATION 
SYSTEMS  RESEARCH  AND  DEVELOPMENT  SERVICE 
SPECTRUM  MANAGEMENT  STAFF 


STATEMENT  OF  MISSION 

The  mniion  of  the  Spectrum  Management  Staff  is  to  assist  the  Department  of  State, 
Office  of  Telecommunications  Policy,  and  the  Federal  Communications  Commission  in 
assuring  the  FAA's  and  the  nation's  aviation  interests  with  sufficient  protected 
electromagnetic  telecommunications  resources  throughout  the  world  to  provide  for  the  safe 
conduct  of  aeronautical  flight  by  fostering  effective  and  efficient  use  of  a natural 
resource-the  electromagnetic  radio-frequency  spectrum. 

This  objective  is  achieved  through  the  following  services: 

• Planning  and  defending  the  acquisition  and  retention  of  sufficient  radio-frequency 
spectrum  to  support  the  aeronautical  interests  of  the  nation,  at  honw  and  abroad,  and 
spectrum  standardization  for  the  world's  aviation  community. 

• Providing  research,  analysis,  engineering,  and  evaluation  in  the  development  of 
spectrum  related  policy,  planning,  standards,  aiteria,  measurement  equipment,  and 
measurement  techniques. 


• Conducting  electromagnetic  compatibility  analyses  to  determine  intra/inter-system 
viability  and  design  parameters,  to  assure  certification  of  adequate  spectrum  to  support 
system  operational  use  and  projected  growth  patterns,  to  defend  the  aeronautical 
services  spectrum  from  encroachment  by  others,  and  to  provide  for  the  efficient  use  of 
the  aeronautical  spectrum. 

•Developing  automated  frequency-selection  computer  programs/routines  to  provide 
frequency  planning,  frequency  assignment,  and  spectrum  analysis  capabilities  in  the 
spectrum  supporting  the  National  Airspace  System.  • 


• Providing  spectrum  management  consultation,  assistance,  and  guidance  to  all  aviation 
interests,  users,  and  providers  of  equipment  and  services,  both  national  and 
international. 
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The  Federal  Aviation  Administration  (FAA)  requested  the  Electro- 
magnetic Compatibility  Analysis  Center  (ECAC)  to  investigate  the 
. effects  of  terrain  on  air-ground  communications  between  air  traffic 

controllers  and  aircraft.  ECAC  produced  terrain-dependent  coverage 
overlays  at  various  altitudes  for  each  Remote  Communications  Air/ 
Ground  (RCAG)  and  Backup  Emergency  Communications  (BUEC)  site  in  the 
contiguous  United  States. 

The  ECAC  Target  Acquisition  Model  (TAM)  was  used  to  produce 
acquisition  overlays  for  each  RCAG/BUEC  site.  Data  pertinent  to  each 
site  was  acquired  from  ECAC  data  files.  This  data  was  verified  and 
updated  by  the  FAA  before  each  site  was  run  through  the  TAM  program. 
Aircraft  altitudes  above  mean  sea  level  and  topographic  data  for 
360  one-degree  radials  around  the  site  out  to  100  nautical  miles  were 
used  as  inputs  to  the  TAM. 

The  acquisition  overlays  were  photo-reduced  over  standard  World 
Aeronautical  Charts.  The  photo-reduced  overlays  were  then  combined 
into  nine  books,  one  for  each  FAA  region  in  the  contiguous  United 
States . 


EXECUTIVE  SUMMARY 
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SECTION  1 
INTRODUCTION 

BACKGROUND 

The  Federal  Aviation  Administration  (FAA)  wanted  to  determine 
how  air-ground  communication  between  air  traffic  controllers  and  air- 
craft was  affected  by  the  terrain  around  the  ground-based  transmitter 
site.  The  Electromagnetic  Compatibility  Analysis  Center  (ECAC)  was 
requested  by  the  FAA  to  determine  the  terrain  coverage  for  all  Remote 
Communications  Air/Ground  (RCAG)  and  Backup  Emergency  Communications 
(BUEC)  sites  existing  in  the  contiguous  United  States. 

OBJECTIVE 

The  objective  of  this  analysis  was  to  produce  contours  that 
depict  where  radio  communications  with  aircraft  can  be  maintained  at 
various  altitudes  around  the  RCAG/BUEC  sites  within  the  contiguous 
United  States. 

APPROACH 


The  ECAC  Target  Acquisition  Model  (TAM),  in  conjunction  with  the 
CALCOMP  plotter,  was  used  to  produce  coverage  overlays.  The  coverage 
included  a line-of-site  contour  and  acquisition  contours  for  aircraft 
altitudes  beginning  with  the  first  even  thousand  foot  contour  above 
the  site  elevation  and  every  2,000  feet  thereafter  up  to  a maximum  of 
20,000  feet.  Terrain  data  was  obtained  from  ECAC’s  topographic  files. 
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SECTION  2 
MODEL  DESCRIPTION 

TAM  DESCRIPTION 

The  Target  Acquisition  Model  (TAM)  was  developed  to  generate 
acquisition  contours  around  specified  site  locations.  These  contours 
are  used  to  assess  site  suitability  or  vulnerability  for  the  detec- 
tion of  aircraft. 

TAM  uses  discrete  terrain  data  to  develop  radar  acquisition  con- 
tour information.  For  a given  site  location,  the  definition  of  the 
terrain  surrounding  the  site,  and  a target  altitude  (above  sea  level 
or  above  terrain),  the  program  generates  terrain  elevation  profiles 
along  a series  of  equally  spaced  radials  emanating  from  the  site. 

The  elevations  are  corrected  for  the  effects  of  atmospheric  refract iv- 
ity  and  the  earth's  curvature.  The  ranges  at  which  targets,  approach- 
ing the  site  at  a given  altitude,  would  first  be  detected  by  radar  is 
determined  by  analyzing  these  profiles.  Ranges  from  each  profile  are 
then  consolidated  to  produce  an  acquisition  contour  around  the  site 
for  the  given  altitude.  APPENDIXES  A and  B describe  the  analytical 
procedures  used  in  TAM  for  profiling  and  line-of-sight  calculations. 

When  scaled  and  reproduced  on  a transparent  plastic  sheet,  the 
acquisition  contours  around  the  site  become  an  actual  map  overlay. 

When  the  overlay  is  properly  oriented  on  an  identically  scaled 
map,  the  contour  lines  for  the  site  indicate  the  distance  to  the 
points  where  an  aircraft  at  the  given  altitude  will  be  within  line- 
of-sight.  FIGURE  1 is  a sample  overlay. 
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SECTION  3 
ANALYSIS 

The  TAM  program  was  used  to  generate  coverage  contours  for  each 
FAA  RCAG/BUEC  site  in  the  contiguous  United  States.  The  latitude, 
longitude,  facility  identifier,  antenna  height,  and  site  elevation 
were  required  for  each  site.  This  information  was  obtained  from  the 
ECAC  data  files.  The  data  was  tabulated  and  sent  to  the  FAA  for 
validation  and  updating  before  being  processed  by  the  TAM  program. 

In  generating  the  coverage  overlays,  the  TAM  program  used  topo- 
graphic data  composed  of  terrain  heights  at  30-second  intervals  in 
latitude  and  longitude.  Inputs  to  TAM  used  to  produce  the  overlays 
were  as  follows: 

1.  Aircraft  altitudes  above  MSL. 

2.  Three  hundred  and  sixty  one-degree  radials  around  the 

site. 

3.  Topographical  data  was  interpolated  every  15  seconds 
along  each  radial  to  a maximum  of  100  nautical  miles. 

The  overlays  were  generated  at  the  same  scale  (1,000,000:1)  as 
the  standard  World  Aeronautical  Chart  (WAC) . The  overlays  were 
photographed  over  the  WAC's  and  reduced.  The  reduced  photographs 
were  arranged  into  books,  one  for  each  FAA  region. 

The  book  for  each  region  was  organized  alphabetically  beginning 
with  the  Air  Route  Traffic  Control  Centers  (ARTCC)  in  the  region. 

The  RCAG/BUEC  sites  serving  each  ARTCC  were  then  alphabetized. 
RCAG/BUEC  sites  within  a region  but  serving  other  region's  ARTCC 's 
completed  the  listing  and  were  organized  in  the  same  manner. 
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In  interpreting  the  coverage  contours,  the  user  should  con- 
sider the  following  limitations.  The  topography  of  CONUS  was 
digitized  at  30-second  intervals.  When  an  elevation  along  a 
radial  profile  lay  between  digitized  topographic  data,  a two-way 
linear  interpolation  was  computed  with  the  four  adjacent  data 
points  to  recover  the  elevation  (see  APPENDIX  A) . In  rugged 
terrain,  the  digitized  data  and  the  interpolation  process  can 
overlook  mountains  or  valleys  that  lay  between  recorded  topo- 
graphic elevations.  Also,  the  topographic  data  file  does  not 
reflect  man-made  obstructions. 

A single  atmospheric  refractivity  index  (310)  was  used  for 
all  of  CONUS.  The  refractivity  index  was  used  to  alter  the  earth's 
curvature  by  increasing  the  earth's  radius.  This  change  of  radius 
compensates  for  the  bending  of  radio  waves  propagated  through  the 
atmosphere.  Along  the  profiles,  elevations  are  also  corrected 
using  the  adjusted  radius  (see  APPENDIX  B) . 

Frequency,  power,  antenna  gain,  etc.  are  not  used  in  developing 
the  TAM  contours.  The  contours  are  based  on  geometric  considerations. 
Consequently,  the  contours  represent  the  maximum  possible  range  a 
target  could  be  detected  at  a given  altitude. 
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SECTION  4 
RESULTS 

RESULTS 

Nine  books  containing  coverage  contours  for  RCAG/BUEC  sites  in 
each  of  the  nine  FAA  regions  in  the  contiguous  United  States,  were 
delivered  under  separate  cover.  TABLE  1 shows  the  number  of  RCAG 
sites  and  the  number  of  BUEC  sites  analyzed  for  each  of  the  nine  FAA 
regions . 


TABLE  1 

NUMBER  OF  RCAG  AND  BUEC  SITES  ANALYZED 
PER  REGION 


Sites  Serving  ARTCC's 
Within  the  Region 

Sites  Within  the  Region 
Serving  Other  Region's  ARTCC's 

Region 

RCAG 

BUEC 

RCAG 

BUEC 

Total 

NE 

4 

1 

0 

23 

EA 

17 

17 

3 

76 

SO 

71 

2 

14 

1 

88 

GL 

120 

8 

6 

0 

134 

CE 

29 

3 

19 

2 

53 

SW 

75 

6 

9 

0 

90 

RM 

55 

7 

13 

1 

76 

NW 

16 

6 

7 

0 

29 

WE 

30 

15 

10 

0 

55 

Notes:  1.  RCAG/BUEC  sites  were  counted  as  RCAG  sites. 

2.  Duplicate  sites,  in  a given  region,  were  only  counted 
once. 
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APPENDIX  A 
PROFILE  GENERATION 


A profile  is  generated  between  two  points  by  calculating  elevation 
values  through  interpolation  of  topographic  data  at  a given  increoent 
along  a great  circle  path  from  the  first  point  to  the  second.  A 
thirty  second  interval  results  in  elevations  being  calculated  approxi- 
■ately  every  half  mile  along  the  path,  while  a fifteen  second  interval 
yields  elevations  at  approximately  every  quarter  mile. 

GREAT  CIRCLE  PATH  COMPUTATION 

Calculation  of  coordinates  along  a great  circle  path  is 
based  on  a unit  sphere  with  the  coordinate  system  shown  in 
Figure  A -1.  The  i,  j,  and  k vectors  are  of  unit  length  in  the 

X,  Y,  and  Z directions,  respectively. 

Transformations  involved  are: 

a.  Spherical  to  rectangular: 

X ■ Cos  ^ Cos  A 

Y ■ Sin  gi  Cos  X (A-l) 

Z ■ Sin  X 

b.  Rectangular  to  spherical: 


X « Tan”^ 
0 = Tan”^ 


CA-2) 


where 

X is  latitude  (positive  north) 
0 is  longitude  (positive  west) 
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Olv«n  a alta  shown  in  Figure  A-2:  as  poin^  S on  the  unit 
sphere  (latitude  longitude  0,),  a vector  S from  the  center 
of  the  sphere  to  point  S is  constructed  as: 


S ■ Sj^  i + Sj  j + Sj  k 


(A-3) 


where  S]^,  M20  and  S3  are  scalars  calculated  from  equation  (a-1) 
with  (Xg,  0a) • Similarly*  a vector  G from  the  center  of  the 
sphere  to  a grid  point  6(Xq,  is  constructed: 

S - 9ri  i + ^2  i + ^3  ^ (A-4) 

The  vector  P from  S to  G is  then 

P ■ G - S - (gj^  - Sj^)  i + (gj  - Sj)  j + Cgj  - s^)  k 

(A-5) 

■ Pi  i ♦ P2  J ♦ P3  X 

The  unit  vector  U in  the  direction  of  P is 


- _ - Pi  i+  P2  3 + P3 

u “ TTT  • = 

VPl^  + P2^  + P3^ 


(A-6) 


* Uj^  i + U2  j + 03  k. 

Equal  spacing  along  the  great  circle  path  from  S to  G is 

achieved  through  the  calculation  of  the  angle  y (shown  in 

n 

Figure  A- 3) as  follows: 

S ■ Cos*^  |sin  Xg  Sin  Xq  + cos  Xg  Cos  Xq  cos  (0q-  0^ 


•^o  - C A, 


(A- 7) 
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ia  user  specified  spacing  in  seconds « 


C is  seconds  to  radians  conversion  factor*  and  other 
variables  are  as  previously  defined. 

coordinates  (Xn*  equal  spacing  along  the  great 

circle  path  from  S to  G are  then  calculated  as  follows 
(see  Figure  A-4)  : 

h * Cos  (^)  (A-8) 

b a Sin  (-j) 
m ■ h tan  (7  - Yn^ 

An  “ b - 

"Cfi  ■ ? + ■ ®ln^  ®2n3  + ®3n^'  1»...*N 

Equation  CA-2)  with  cin*  C2n»  C3n  then  gives  Xn/  ffn. 

These  coordinates  apply  equally  to  the_  unit  vector 
(Figure  A-3),  in  the  direction  of  C^,  which  is  the 
vector  from  the  center  of  the  unit  sphere  to  the  nth 
point  on  the  great  circle  path  from  S to  G. 


DATA  RETRIEVAL 

It  is,  of  course,  necessary  to  determine  which  topographic  data 
block  is  needed  to  obtain  the  elevation  of  a given  point,  and  then 
to  retrieve  it  from  drum.  In  order  to  minimize  drum  access  time, 
two  data  blocks  are  always  retained  in  computer  memory  during  profile 
generation  from  a given  site.  One  array  is  reserved  for  the  block 
area  containing  the  site  location.  The  second  array  is  used  to  store 
other  data  blocks  which  are  required  for  grid  points  not  contained  in 
the  site  block  area  (See  FIGURE  A-5) . 
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Let  X ifi  be  the  latitude  and  longitude  (in  radians)  of  the  n 
point  alonj  great  circle  path  from  the  site  to  a given  grid  point. 
Since  topographic  data  blocks  are  referenced  by  latitude  and  longitude 
of  the  southirest  comer  in  ainutes,  the  block  reference  (Xj^,  for 

X , <fi  and  additional  measures  (x,  y)  used  in  elevation  extraction  are 
determined  as  follows: 


For  >0, 

^ 

y • Xg  - 

For  Xjj  < 0, 

X«w  • - I*V 

y * Xg  - + Ly 

For  «pn  > 

<P*w  ■ ‘P'T  ^ 

X • <pt  “ ^Pc 

For  ©n  * 0, 

'Psw  * ‘PT 

X • ©Y  - ©g 

whore 


(A-9) 


^sw*  latitude  and  longitude  of  the  southwest 

comer  in  minutes, 

Xc/  ©c  latitude  and  longitude  converted  to 
minutes , 
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. Lyin  minutes. 


. Lgin  minutes. 


Ly,Lg  are  vertical  and  horizontal  topographic  block  lengths 
in  minutes, 

X,  y are  extraction  measures  in  minutes  (see  Figure  A-6), 
and  (>  denotes  truncation  to  nearest  integer  value. 


The  necessary  block  reference  point,  (Xsw*  <Psw) « is  then 
compared  wi^  those  of  the  two  data  blocks  in  the  computer,  if 
the  required  data  is  already  available,  the  elevation  at 
On  Biay  then  be  extracted;  if  not  present,  the  data  is  first 
retrieved  from  drum  storage. 

ELEVATIOM  EXTRACTION 

Two  methods  of  deriving  an  elevation  are  available.  The 
first  procedure  assumes  the  same  elevation  for  the  location 
in  question  as  is  recorded  in  the  topographic  data  file  for  the 
location  nearest  to  it.  The  second  method  performs  a linear 
interpolation  in  two  directions  using  the  elevations  of  the 
four  proximal  points,  with  the  resultant  being  the  desired 
elevation. 

The  required  row(s)  and  colvunn(a)  of  the  121  x 121  topo- 
graphic data  matrix  (see  Reference  1)  are  first  determined. 


m - <M> 
n - <N> 


(A- 10) 


where 


m is  the  number  of  the  nearest  matrix  row  south  of  the 
location  in  question, 

n is  the  number  of  the  nearest  matrix  column  west  of  the 
location  in  question. 
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N = 1. 

hi 


and  other  variables  are  as  previously  defined. 

The  relation  of  the  proximal  elevation  shown  in  Figure  A-7  to 
the  desired  elevation,  E^,  is  determined  by  the  following  weights: 

Aj  = N - n 

A2  = 1 - Aj 

6^  = M - m 

62  = 1 - Aj 

The  desired  elevation,  E^,  is  then  extracted  by  one  of  two 
methods : 

1.  Nearest  Point  Option. 

If  6^  < 62,  i = m;  otherwise,  i = m + 1; 

If  Aj  < A2,  j = n;  otherwise,  j = n + 1. 

The  desired  elevation  is  then  obtained. 


E,  = H.  . 
d i,j 


(A-11) 


2.  Four  Point  Option. 


Interpolation  in  the  north-south  direction  gives 

^a  ^ *^2  \,n  '*1  ^m+l,n 

^b  ^ "^2  \,n+l  “^l  ^m+l,n+l 


(A- 12) 


Interpolation  in  the  east-west  direction  then  results  in  the 
desired  elevation. 


^d  ^2  * ^1  ^b 


(A- 13) 
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Since  elevation  values  are  initially  biased  by  a quantity  AE, 
quantized  by  a factor  Q and  stored  in  compressed  form  in  the  topo- 
graphic data  file,  the  elevation  derived  from  either  method  must 
be  converted  to  feet  above  sea  level  by  the  equation: 

E = Q + AE  (A- 14) 

The  extracted  elevation,  E,  is  thus  one  of  any  number  of  elevations 
that  constitute  an  elevation  profile  from  a site  to  a given  grid 
point . 


L 


20 


FAA-RD-77-126 


Appendix  B 


APPENDIX  B 


LINE-OF-SIGHT  DETERMINATION 


B®Fore  line-of-sight  c&n  be  detenoined  along  a profile,  two 
factors  oust  be  taken  into  consideration:  the  effect  of  refractivity, 
or  bending,  of  radio  waves  as  they  are  transmitted  through  the  at- 
mosphere, and  the  effect  of  earth's  curvature  with  respect  to  the 
elevations . 

Atmospheric  refractivity  is  accounted  for  by  assuming  an  ef- 
fective earth's  radius,  a = ka  where  k is  the  effective  earth's 
radius  factor  and  a is  thi  true  radius  of  the  earth.  This  allows 
radio  waves  to  be  drawn  as  straight  lines  over  an  earth  with  ef- 
fective radius  a , 
e 

The  value  of  k,  for  rays  leaving  a transmitter  at  low  angles 
to  the  horizon,  is  given  by  the  expression 


where 


1 * IQ-^N 
1 + 10"^N  + 10‘® 


local  refractivity  index 
true  earth's  radius 


(B-1) 


AN  = change  in  refractivity  with  respect  to  altitude. 

In  general,  the  change  of  refractivity  varies  exponentially  and 
is  given  as  AN  = -7.32  e0-005577N,  Taking  a = 6370  kilometers  yields 


1 * IQ-^N 

1 + 10‘®  N - 0.04663 


(B-2) 


^suming  that  the  effective  earth's  radius  is  ka,  each  elevation,  En, 
is  now  reduced  by  an  appropriate  amount,  Ahn,  which  corrects  the 
elevation  for  earth's  curvature.  Using  Figure  B-1,  it  is  shown  that 


Ahn  = -ka  + ^(ka)^  + Xn^ 


(B-3) 
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The  arc  length,  dn,  is  a reasonable  approximation  for  Xn  as  Xn 
differs .from  dn  by  less  than  one  mile  for  path  lengths  up  to 
425  miles.  A binomial  series  ejqpansion  of  the  radical  in  Equation 
B-3  yields  the  following  equation 

Ahn  - (B-4) 


The  corrected  elevation,  Ec,  is  then 

_ dn^ 

Ec  » en  - ^ 


(B-5) 


The  vertical  angle,  0n,  is  then  calculated  (Figure  B-2)  from 
the  site  to  each  corrected  elevation  along  the  profile  by  the 
following  equation: 


en  = 


Ec  - CEs  * H) 
dn 


(B-6) 


where 


Es  = site  elevation  above  sea  level  (feet) 

H = antenna  height  (feet) . 

As  shown  in  Figure  B-3  the  limiting  line-of-sight  will  occur  at 
the  values  of  dn  and  Ec  which  correspond  to  the  maximum  value  of  the 
vertical  angle  9n,  or  the  radio  horizon  angle. 


FIGURE  B-1.  ELEVATION  CORRECTED  FOR  ATMOSPHERIC 
REFRACTIVITY  AND  EARTH'S  CURVATURE. 
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The  Line-of-Sight  determination  described  in  this  appendix 
assumes  that  radio  wave  bending  has  a linear  dependence  on  the 
refractive  index.  For  heights  above  one  kilometer,  the  exponential 
dependence  on  height  should  be  e^iloyed  in  lieu  of  the  linear  ap- 
proximation. This  is  acco^>lished  by  using  the  NBS  Ray  Trace 
Correction. ^ 


^Bean,  B.  R.,  Thayer,  G.  D.,  CRPL  Exponential  Reference  Atmosphere, 
NBS  Monograph  4,  October  29,  1959. 
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